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Abstract: 

The time evolution of the coherent forward scattering of the synchrotron 
radiation for resonant nuclei in Brownian motion is studied. Apart from 
target thickness, the appearance of the dynamical beats also depends on 
'a' which is the ratio of the harmonic force constant to the damping force 
constant of harmonic oscillator undergoing Brownian motion. 
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Introduction: 

In recent years , synchrotron radiation (SR) has been used to study 
Mossbauer spectroscopy in the time domain. By measuring the time de- 
pendence of the intensity of nuclear de-excitation [1 and references therein] 
various solid state and nuclear parameters are being studied. Synchrotron 
radiation is produced when ultra-relativistic electrons or positrons acceler- 
ate while transgressing bending magnets or periodic magnetic structures and 
emit intense radiation. SR is emitted in short repetitive pulses whereas emis- 
sion of photons from a radioactive nucleus are random in time. 

The Mossbauer spectra of iron-protein crystals have special shape which 
provide a unique signature, indicating that Mossbauer atom undergoes 
bounded diffusive motions in the biological systems. It is observed for biolog- 
ical systems that above a particular temperature, which is the characteristic 
of the system (i)the Mossbauer lineshape becomes non-Lorentzian, (ii) the 
mean square displacement (msd) shows temperature dependence which is not 
typical of either the Einstein or the Debye model and (iii) the width of the 
Mossbauer lineshape increases with the temperature. The non-Lorentzian 
lineshape was observed by Parak et. al. [2], Cohen et. al.[3] in myoglobin. 
Similar phenomena was seen in haemoglobin by Mayo, Parak and Mossbauer 
[4]. Pioneering studies on the dynamics of proteins using Mossbauer Spec- 
troscopy were done by Keller and Debrunner [5], Parak et. al. [6,7]. These 
unusual results observed for biological systems, have been explained by us- 
ing a dynamical model of the bounded diffusion. These models are based 
on the concept of Brownian motion of an overdamped harmonic oscillator 
[2,3,8,9]. Recently synchrotron radiation has been used to study dynamics 
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of Mossbauer atom in the biological systems [10,11,12,13]. The goal of the 
present paper is to study the time evolution of resonant scattering of SR from 
nuclei undergoing Brownian motion in harmonic potential. 
Methodology: 

In this study we use mathematical treatment developed in reference [8,9] 
as a starting point. A short pulse of SR can be decomposed into a contin- 
uum set of monochromatic spectral components within a frequency interval. 
Smirnov and Kohn [14,15,20] calculated the forward transmitted wave packet 
by integrating over all spectral components of the transmitted radiation to 
obtain 



E(t,Z) is the time dependent electric field amplitude of the synchrotron 
radiation transmitted through a nuclear target of thickness Z. The function 
E u0 (Z) has the modulus (^j) 1 ^ 2 exp(— ^) . The electronic absorption 
coefficient is given as [i e =K \ where x = ^Ri °e being the total cross section 
of electronic inelastic scattering and Vq is the target volume corresponding 
to one nucleus. K is the wave number defined as K=^ and uj is the frequency 
having wave vector k. I is the intensity of the SR within the frequency band 
Auj, as determined by the monochromator system. C(k, uj) is the nuclear part 
of the susceptibility of the target and is related to the scattering amplitude 
and can be written as 
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characterizes the strength of the nuclear response at the resonant frequency 
uj eg and (f)(k,cu — u> eg ) is the universal resonance function. Here,g and e define 
sublevels of the ground and the excited states of nucleus respectively and To 
is the natural linewidth of the excited levels, /lm(^) is the Lamb Mossbauer 
factor, I g is the nuclear spin in the ground state while < g\j (k)\e > rep- 
resents the matrix element of the scaler component of the nuclear current 
density vector. 

At resonance , nuclear part of the susceptibility of the target has a fre- 
quency dependence , which is determined by the the nature of universal 
resonance function, obtained (in the present case) by taking into account the 
diffusive motion. It is described by the following 

<j)(k,uj)=J o dtexp(-iut--£r)F a (k,t) (4) 

where 

F s (k,t) = J drexp(-ik.r)G(r,t) (5) 

Here G(r,t) is the correlation function. The derivation of above equations 
is given in references [14,15] in detail. We consider only that case where 
polarization is absent in the coherent forward scattering. The most frequent 
application of SR in the nuclear resonance spectroscopy is the measurement 
of the time-dependence of the forward scattering intensity, given by 

I f (t,z) = \E(t,Z)\ 2 (6) 

It is interesting to note that Mossbauer absorption spectrum is dependent 
on the real part of the universal resonance function <p(k,uj). On the other 
hand, the time dependence is determined by the entire universal function. 
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Harmonically bound nuclei in Brownian motion: 

The time evolution of nuclear resonant scattering of SR pulse by nuclei 
undergoing various types of diffusion (like free diffusion, continuous localised 
and jump diffusion) has been studied earlier [14,15] . Here we consider an 
important case of diffusion in the form of Brownian motion. Gunther et. al. 
[16,17] were the first to predict that Brownian motion of a Mossbauer atom 
results in the non-Lorentzian nature of Mossbauer line shape. The theory for 
the cases of motion of harmonically bound Mossbauer atom in the Brownian 
motion has been extensively developed both for the classical [2,3,8,9] and 
quantum cases [18]. In the present paper we will confine to the classical case 
and use the theory as developed by Nowik et. al. [8,9]. For one- dimensional 
harmonic oscillator in the Brownian motion, Uhlenbeck and Ornstein [19] 
derived a general formula for self correlation function G(x, xq, t) which is the 
probability that at time t the nucleus will be at position x if at time t=0 it 
was at position x . This self correlation is given as [8] 

a , 1/2 ,—a(x - x exp(-at)) 2 , 

G(x,x ,t) = ( ; ; -) 1/2 exp( K - — P (7) 

v ' °' ; v 2tt D(l - exp(-2at)' PV 2D(1 - exp(-2at)) ' K ' 



where 



w 2 



* = J (8) 



is the ratio of the harmonic force constant w and damping force constant (3. 
The diffusion constant D is given by 

D=±% (9) 

Nowik et. al.[8,9] generalised this self-correlation function to three dimen- 
sions for the over-damped case and obtained a simple formula for the line 
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shape. Mossbauer spectra were computed for range of a values and fixed dif- 
fusion constant D. The parameter a [8,9] decides the nature of the Mossbauer 
lineshape. The Mossbauer line-shape is non-Lorentzian for smaller values of 
a [8,9]. Since Mossbauer lineshape for biological systems is non-Lorentzian 
in nature , we will focus on smaller values of a in our calculations. 

The universal resonance function characteristic of the Brownian motion 
can be written as [9] 

ftk, U ) = exp(— ) V —A—T-, U ° . r (10) 



where 



Thus, in the present case, for equation (1), C(k,u>) is given by 

C(M) = EE7 B ^ m) r (12) 

where 

a k (m) = —A—) m ezp(-— ) (13) 
ml a a 

Thus E(Z,t) is given as 

E(Z, t) = %^ [ du exp{-iut) exp (- M™) ) (14) 

where 

^ = X£B ge (15) 

The summation over 'ge' in equation (15) is related to the summing over 
the transition between the ground and the excited states of nucleus [ 14,15 
]. It is important to note that for any particular value of (^p) ,only a few 



Lorentzians contribute to the universal line shape [9]. 
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Results: 

To illustrate the role of Brownian motion in the time dependence of for- 
ward scattering,we first consider of a thin target (single line sample) 
having an effective resonance thickness fi n Z=l. The results of the numerical 
calculation for various values of a are shown in the Figure 1. It is evident 
from this figure, that the natural decay appears (manifested by straight seg- 
ment of time response curve) faster for smaller values of a and slower for 
larger values of a. The time response in the case of single-line target, having 
an effective resonant thickness fi n Z=10 (thick samples) for various values of 
a, are shown in the Figure 2 where dynamical beats appear. It is also clear 
from this figure that smaller the value of a, earlier the dynamical beat occurs 
and vice versa. The appearance of the dynamical beats is also dependent on 
the value of [i n Z. Comparing Figures 1 and 2 it is clear that a and fi n Z share 
a definite relationship, with each other so for as the appearance of dynamical 
beats are concerned. In order to further investigate this relationship we fix 
the value of a and see how the variation of fi n Z affects the appearance of 
dynamical beats in the time response behaviour. The results are shown in 
the Figure 3 which indicates that larger the value of fi n Z, earlier does the 
dynamical beat appear and vice versa. 
Comparison with diffusion results: 

In general the Mossbauer lineshape has a Lorentzian character. The diffu- 
sion of the Mossbauer atom causes broadening in this Lorentzian lineshape. 
For the case of free diffusion, the increase in the Mossbauer linewidth is 
proportional to the diffusion constant. The presence of the this diffusion 
constant in the exponential factor of the time response causes an accelerated 
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decay of the coherent signal. However, the position of the appearance of 
the 'dynamical beat minimum ' is independent of the value of the diffusion 
constant. All dynamical beats appear at the same time with variable mini- 
mum values for different values of the diffusion constants. But in the cases 
of bounded diffusion inside a potential well and jump diffusion, the universal 
resonance functions have complicated shapes represented, in general, by the 
coherent superposition of the Lorentzian lineshapes. The width and weight 
of each Lorentzian is determined by the specific nature of the diffusion pro- 
cess. Thus, the appearance of the dynamical beats in the time response is 
sensitive to the diffusion rate and cage-size for the cases of the jump diffusion 
and bounded diffusion,respectively. 

The universal resonance function for the case of Brownian motion in gen- 
eral (and for smaller values of a in particular) is non-Lorentzian.This is be- 
cause the universal resonance function is the result of superposition of vari- 
ous Lorentzian line shapes. Each Lorentzian line shape has its characteristic 
weight and width as determined by the value of a . This type of the nature 
of universal resonance function results in the complicated behaviour in the 
time response. As in the cases of continuous diffusion and jump diffusion, 
the physical reason for the complicated behaviour in the time response for 
the case of Brownian motion,is the split of the universal resonance function 
into several terms. 
Conclusion: 

The method of coherent forward scattering of the SR reveals complex 
amplitudes of the oscillation of the electromagnetic field in the presence of 
the Brownian motion. 
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Figure captions: 

Figure 1: The time dependence of nuclear forward scattering of the syn- 
chrotron radiation in the presence of Brownian motion for different values 
for a thin sample of effective thickness [i n Z=l corresponding to a = 0.5, 0.6, 
0.7, 0.8, 0.9, 1.0, 2.0, 4.0 . The lowermost curve corresponds to a=0.5 and 
uppermost curve corresponds to a=4.0. a is units of sec -1 . 

Figure 2: The time dependence of nuclear forward scattering of the syn- 
chrotron radiation in the presence of Brownian motion for different values of 
a for a thick sample of effective thickness [i n Z=10. a is units of sec -1 . (1) 
a = 0.5;(2) a=0.6;(3) a=0.7; (4) a=0.8; (5) a =0.9; (6) a =1.0; (7) a =2.0 
; (8) a=4.0 . 

Figure 3: The time dependence of nuclear forward scattering of the of 
synchrotron radiation with a=l for different values of fi n Z. a is in units of 
sec" 1 (1) /i n Z = 50.0;(2) /i n Z=40.0;(3) /i n Z=30.0;(4) /j n Z=20.0; (5) fi n Z 
= 10.0; (6) fi n Z =5.0 ; (7) fi n Z =1.0 . 
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